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I, Introduction 

The classical root locus has proven to be a valued^Xe analysis and 
design tool for single input single output linear control systeas. 
Heseeuxh is currently underway to extend these methods to laultt-input 
Qulti^output linear control systems and linear optizaal control systems* 

In this paper wc present analysis techniques for both of these multi- 
variable root loci. We show how to compute angles and sensitivities for 
both types of loci, and how to determine the asymptotic behavior (as 
control weights get small) of the optimal root locus. 

Previous work on angles and sensitivities is contained in [1#2] • The 
former uses time domain techniques (the eigenvalue problem) and the 
latter uses frequency domain techniques. We extend the time domain 
techniques through the use of generalized eigenvalue problems, and 
wo show this approach to be significantly better for computing angles 
of approach. 

Perhaps the most significant development in the understanding of 
multivariable root loci was the concept of multivariable transaission 
zeroes [3] • These form the endpoints of all asymptotically finite 
branches. Determining the behavior of the asymptotically infinite 
branches, however, has proven to be a difficult problem and all of the 
details are not yet Icnown [4] . Frequency domain interpetations of 
multivariable root loci using Riemann surfaces have been given (5], 
and the behavior of the closed loop eigenvectors has also attracted 


some attention. 
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The root loci of linear quadratic optijaal control syetena were 
first described for single-input single-output systeas in [6,7]. These 
methods have been extended to the nmlti-lnj^t case in [8,9,103. Asymptotic 
properties (which include the asymptotic behavior of eigenvectors) are 
used to select quadratic weights In (llj . Optimal root loci can be 
considered a special case of ordinary linear feedback loci, and it 
turns out that the asymptotically infinite behavior of this special 
case is better behaved. Consequently more progress has been made in 
analyzing this behavior [12,13J. He extend the available analysis 
techniques for determining the asymptotically infinite behavior to 
Include the behavior of the eigenvectors. In doing so we use a new 

type of subspace decomposition which simplifies the previous analysis 
technique (12J . 

In section II we develop the formulas for computing angles and 
sensitivities of the multivariable root locus. In section III these 
formulas are applied to the multivariable optimal root locus. Then in 
section IV we develop analysis techniques for deteraing the asymptotically 
infinite behavior of the multivariable optimal root locus. 

The Generalized Eigenvalue Problem 

The generalized eigenvalue problem is to find all finite X and 
their associated eigenvectors v which satisfy 

Lv " iMv 

I. and H are real valued p>:p matrices which are not necessarily full 
rank, if M is invertible then prcmultiplication by changes the 
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generallzea eigenvalue problea into a standard eigenvalue problea for 
»rtiich there are exactly p solutions. Ii. general there are 0 to p finite 
solutions, except for the degena’-ate case trften all A la the complex plane 
is a solution. Reliable FORTRAN subroutines based on stable numerical 
algorittens exist in EISPACK [141 to solve the generalized eigenv^auo 
problem, see [151 for the application of this software to a related 
class of problems. 

Notation 

Hatrices are denoted by capital letters, scalars and vectors by lower 

case letters. and y” are the transpose and Hermitian transpose, 

respectively, of A and y. a‘*^ indicates or, equivalent, (a’^)’^. 

A > 0 and A > 0 indicates that A is positive seraidefinite and positive 

definite. If A is synnetric then A^''^ is the (nonunique) decomposition 
1/2 1/2T 

of A into A A . iTubspaces sire denoted by script letters, with 
the exception of the real vector space a”, "im A" and "ker A" are the 
image and kernel of the linear nap A. The dimension of U is dim O, 
subspace Inclusion is C subspace intersection is O, and a linear 

combination of subspaces in U + V. An open loop linear system is denoted 
by (A,B,C). 
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II, Angles and Sensitivities of the Root Locus 

We consider the linear tiae invuiant output feedback problea 


X > Ax •<> Bu 

a IS 

X e 21 , u C 21 

(1) 

y “ Cx 

a 

y CB 

(2) 

o-^rsr . 


(3) 


The closed loop systea matrix and its eigenvalues, right eigenvectors, and 
left eigenvectors are defined in the usual way by 

(4) 

° i “ 1,..., n (5) 

s^I) - 0 i » 1,..., n . (6) 

We make the assumptions that (A,B) is controllable, (C,A) is observable, 
and K is invertable. Only the case where the number of inputs and outputs 
are equal is treated, and we further assume that the closed loop eigen- 
values axe distinct. 

As k is varied from Infinity down to zero the closed loop eigenvalues 
trace out a root locus. At k - * (to be more precise let I » lA awi use 
it *» 0) the n branches of the root locus start at the open loop eigenvalues. 

k 0, some number p < n-m of these branches approach transmission 
zeros, which are defined here to be those values of s which reduce the 
rank of 

A - si B 

- C 0 
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Wo further rule out degenerate cases ^ in other words we assume that 
and C do not conspire in such a way that any value of s in the complex 
plane reduces the rank of the matrix. V7e rate that if all of our assumptions 
are valid then this is an adequate definition of transmission zeros [ 3 ] • 

Also as k Of the remaining n-p branches of the root locus approach in- 
finity. The behavior of these branches concern us in Section r/. 

At any point on tlie root locus an angle can be defined. Consider the 
closed loop eigenvalue which is computed for some value cf k. If k 
is perturbed by an amount Ak then will be perturbed by A»^. As 
Ak 0 then As^/Ak approaches the constant ds^/dk (if this limit exists) . 

The 2uigle of the root locus at poinr is then defined to be 

^ - arg(ds^)f 


^ere ”arg" is the £Lrgument of a complex number. The angles of the root 
locus at the open loop eigenvalues are the angles of departuref and the 
angles at the transmissions zeroes are the angles of arrival^ rigure 1 
illustrates these definitions. 

Next we define the sensitivity of a closed loop eigenvalue to a change 
in k to be 


This definition is motivated by the approximation 




f 


from which we obtain 




I X OPEN LOOP POLES 


O TRANSMISSION ZEROS 



j-PLANE 


OF DEPARTURE 
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So, to first order, a change Ak will laove a^ a distance jdo^j In the 
direction arg (da^) . 

Before presenting formulaa for these angles and sensitivities, we 
present the following lesm, which shows how the generalized eigenvalue 
probleo can be used to compute the closed loop eigenstruoture. 


Sj, , x^, and are solutions of the generalized eigen* 

value probleas 



Btoof , Prea (7) we see that 


(A - + BWj. ■ 0 




Substitute (10) into (9) to get 


(A - s^Dx^ - ^ BKC » 0 , 


( 9 ) 

( 10 ) 


which is the same as (5), the defining equation for the closed loop eigen- 
values and right eigenvectors. In a similar way (8) can bo reduced to (6) 
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I 


) 


and the proof is coaplete. Letnaa 1 Is not a new result but we have been 
unable to find a reference for it, A precursor of this result (without 
consideration of closal loop eigenvectors# and without isention of the 
generalized eigenvalue problen) is the polynonial systea natrlx repre- 
sentation of Rcsenbrock 13]. 

When k>0 then p, the number of finite solutions of s^ in (7) and 
(8), is equal to n. When k»0 then 0 < p ^ n-m (under stated assumptions). 
The ability to use (7) and (Q) with k= 0 is the major advantage of the 
generalized eigenvalue problen. The finite solutions s^ when k- 0 are 
the transaission zeroes of the systaa# and the and y” vectors are the 
right and left zero directions £15] . From (7) and (8) it is clear that 
as k -► 0 the finite closed loop eigenvalues approach the transmission 
zeroes and the associated eigenvectors approach the zero directions.^ 

The solutions of the generalized eigenvalue problems contain two 
vectors and vdiich do not appear in the solutions of the ordinary 
eigenvalue problems. The importance of the vectors can be explained 
as follows. The closed loop right eigenvector is constrained to lie 
in the m dimensional subspace ofSl" spanned by the columns of (s^l-A)"^B 
117] . Exactly where x^ lies in this subspace is determined by , via 
Xf - (s^I - A) ^BV^. This follows frcra the top part of (7) . If the state 


In (16) , transaission zeroes are computed by solving an eigenvalue problem 
for equation (5) with k close to zero. This is the high gain feedback 
method. In '15] this is shown to have the potential to be computationally 
inferior to solving equation (7) with k =« 0, ^ 


I 
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o£ the closed loop systen at time zero la * « 

g then the state trajectory 

for tl« ,ro.t« 1 , 

u<t) . -(o/klJCjO ‘ . ,Vj.'* . This follow, froo tho botto. part of (7). 
The vectors play an analogous role In the dual system C*, B*^). 

For our purposes, however, tho and nj vectors are also significant 

because they can be used to compute the angles on the root locus. This Is 
shown by the following Tbeoren: 

Theorem 1. The angles of the root locus, for 0 < k < • and for distinct 
are found by 


-y^BRCx, 


arg (ds^) « arg — i— 0 < k ^ » 

i"l,..., p 


Vi 


( 11 ) 




arg (ds_j^) « arg i. q < k < « 

Vi ^ , 

* “ if*ee#P 

R^ark. The angles of departure are found using (11) with k - the 
angles of approach are found using (12) with k- 0. For k > 0, p - n; and 
for k = 0, 0 < p £ n-ffl. 


( 12 ) 


Proof. The proof of (11) can be found in [1] . t 
show that 

ds -n^K'^V 

— i. . — i i 

dk 

Rewrite (7) as 


To prove (12), we first 


Vi 


(13) 


(L - ** 0 


where 


i ** 1#»«*^ p 



"* • 


r ^ 




A B 


I 0 


X . * 

L » 

-1 

t M « 


f V, « 

i 


•C -kK 


0 0 

i 

V 






i 

Let also. 




. • 


( 14 ) 
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< ■ Kn."] 

Differentiate (14) to get 

[dir + (r.-o^H) 0 . 

Multiply on the left by u® to get 

“l[dk - 0 . 

After substituting and rearranging the result is (13). The formula for 
the angle is shown from (13) to be 


arg(ds^) » arg(dk) 


+ arg(-l) + arg j -i- ^ 

/ 


Slnc» k varle, fr„ . t. 0, arg (dk, . igo-, and nine. 

•rg(-l) . leo-, tha „»lt l. (ij, . gni, ceputas kha proof. 

The angles on the root locus for 0 < k < ® can a 

<t s. can be found using 

either (11) or (12) vvfar^*- e < 

w tiaj j a Except for k very close t« j. 

«.jr ciose to zero, when numerical 

problems may be a factor, it is best to use ( 11 ) because it involves 

-Lin, an „ dfoanafonaX afganvaXna prokXe. rakkar kkan an nk. dkaanakonak 
Sdnaralirad alganvalna probl„. gk. 

Obtained from (7) and (8) of Lemma 1, can be used to pass back and forth 
from ( 11 ) to (12): 

Cx^ - -kK”^v^ 

yl’B » kn®Jc-\ 

.ra, kkea. idankikkaa va aaa kkak akan k- 0. Cr, . „ and pja . 0. .kara- 
^ora. ,n, aannok ka „sad wkan k-o ko oon^ra anglaa of approack kao.aaa 
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arg (ds^) “ arg (0), which is not defined.^ 

Equations (11) and (12) are still valid when the controllability 
and observability assumptions are relaxed. However, angles ca n only 
be computed for modes that are both controllable and observable because 
only these modes move as a function of k, and thus have well defined angles. 

n 

The and vectors are also useful for the calculation of eigen- 
value sensitivities* This is shown in the next Iczma* A separate proof 
of this lenaaa is not necessary because the proof follows from intermediate 
steps in the proof of Theorem 1. 


Lcacna 2. The sensitivities of distinct closed loop eigenvalues to 
changes in for 0 £ k ^ are found by 




y^KCx^ 

dk 

.2 

H 


k 

1 


0 < k 

1 « 1, ...rp 


dSi 

n 


dk 


H 





0 £ k < « 

1 ” l^«.*fp 


(15) 


(16) 


Equations (15) and (16) give the same anvers for 0 < k < Even 
though k appears only in (15) , ve note that both (15) and (16) are 
dependent on k because the vectors n”, and are all dependent 

on k* 


In [1] a limiting argument as kr^O is used to derive alternate equations 
for the angles of approach* These results are more complicated than (12) 
because the rank of CB must be determined. The generalized eigenvalue 
problem eliminates the need for this rank detemination. Furthermore, 
the equation given in [1] for the Rank (CB) » m case (3.16b (!]) is 
incorrect due to an error in the derivation after (3.15 [1]). This 
error leads to the incorrect conclusion that the angles of approach 
are independent of the output feedback matrix K. 



r 
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Exarole 1 

st.ate that above resulta, we define a system S(A,B,c) and 

Plot r«t iccl tor „ch of 3 output f.oattOc ottrloeo K. Th, ,pot„ 
niatrices are 


C - 


-4 

0 

4 

0 

C 

8 


7 

3 

7 

-1 

-5 

•14 


•1 

0 

•4 

0 

2 

0 


The outfit feedback matrices 


13 

2 

8 

0 

•2 

2 

are 


0 

1 

2 

•2 


1 

0 

0 

0 


Case # 1 


-r: :i 


Case 


c :i 


Case #3 


K = 


P °] 

LO 50 J 


case ,2 io tho o«o „ „3td i„ 111 . Tho root loci a„ ahoun to 

2. Tho aogloo of Ooporturo ana approach ucro conputod and arc Itotoo 
in Table 1. 

Tho ayoto, haa two open loop uostabla Moa that arc attracted to 
unatabla traoersaioo roro.., ,0 for all values of h the systtr. 1. 
uostable. The systre ha, t„, open loop stable uodes that are attracted 
to - aloe, the negative real one of the branches first goes 

to the right along th. negative real axis and then turns around. 

The turn around point is called a branch point. The rent locus can 
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b. thc,»,ht of OS taio, plotted on . pieoann tnrfnco, „d th. branch 
points are point, at vMch th, root loon, „ov„ batwcan dl££„,„t ah,,t. 
of the Rienann stirface [51. 


table 1 


Angles of Departure and Appr oach for Exanpl» i 




f 


I 


I 

I 

i 


Case 




Case ^2 




ric7»irc 2 



Systi'n with Outinit K.«t\iLMcK 
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XIX* Angles and Scnsitivitios of the Opt imal Root t^ci^ 

Our attention now shifts from the linear output feedback problem to 
the linear optloal state feedback problem with a quadratic cost function. 
As in 17, 121, we show that the optimal root locus for this problea is a 
special case of the ordinary output feedback root locus. We show how 
to compute asimptotically finite properties of the optimal root locus 
and how to compute angles and sensitivities. 

The linear optimal state feedback problem is 


ik «» Ax + Bu 
u ■ F(x) . 


X e K*', 


U 6B 


m 


(17) 

(18) 


The optimal control is required to be a function of the state and to 
minimise the infinite time quadratic cost function 


<10 


pu"^Ru)dt , 


where 

Q “ ^0 

R - > 0 

0 < p < “ 

AS usual we assume that (A,B) is controllable and that the state weighting 
matrix I factored into 


Q - H H , 

whore Rank (Q) - Rank (H) - r, and H e produces an observable 


pair (H.A). 



f 


Kalnan (isi has shown (for p>0) that the optinal control is a linear 


function of the state 


u » -Fx # 


where 


P « — R BP 
P 


and P is the solution of the Riccati equation 


0 - Q + aS + PA - ^ PBR"^B^P 


The closed loop system matrix is 


As p is varied from infinity down to zero the closed loop eigenvalues 
trace out an optimal root locus. 

To study the optimal root locus we define a linear output feedback 
problem with 2n states, m inputs, and m outputs. 


A 0 


-Q -A 


A. T ""I 

C - to D 1 K - R 


The closed loop system matrix is 


r 1 -1 T 

A -iBK B 

Z - A - ^ BKC - „ 

• P _Q _ 
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The Z xnatrix is soaetimes called the Hamiltonian matrix. Its 2n 
eigenvalues arc known to be svmmetric about the imaginary axis, and 
those eigenvalues in the left lialf plane (t«HP) are the eigenvalues of 
We again assume that the 2n eigenvalues of Z are distinct. 
Then the right and left eigenvectors of Z can be defined to be 

n 

i ■ 1, . 2n . 

The right eigenvectors c&n be further deconposed into 



Then the vector is a right eigenvector of and There 

is api«urently not a slailarly easy way to find the left eigenvector 

H H 

y. of A ft from z, and w*. 

■'i cZ i i 

The closed loop eigenvalues, right and left eigenvectors of Z can 
be found by solving ordinary eigenvalues problems. Alternatively, 
using Lemma 1, they can be found by solving the following generalized 
eigenvalue problems. 



i»l,...,2p (23) 



» 0 


i»l,...,2p (24) 


The number of finite generalized eigenvalues is 2p - 2n if p > 0 and 


is 0 < 2p < 2(n-m) if p » 0. 
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if 

Sf 


if 


c 

» 


Wa can analyze the optimal root locus by using the LHP portion 
of the root locus of the Hamiltonian system. At p ■ « the n branches 

of the optimal root locus start at the stable open loop poles (or the 
mirror image about the Imaginary axis of the open loop unstable poles) . 
The branches of the optimal root locus always stay in the LHP. As 
P^O, p of these branches stay finite and approach transmission zeroes, 
idiere 0 ^ p ^ n-m. Those transmission zeroes are the finite LHP solutions 
of the generalized eigenvalue problem (23) with p-0. The right zero 
associated with the transmission zeroes eure the portions 
of the associated z^ vectors.^ 

These asymptotically finite properties will be grouped together 
in the following way: 

S® ■ dlag (Sw ...» s®) 

^ P 

X° - IxJ, ..., x°l . 
i P 

0 0 
Each s^ is a transmission zero and each x^ is a right zero direction. 

Because each x^ is a direction it is only unigue to within a scalar 

multiple. 


If Q « H H, where H S 3R , then n+m dimension generalized eigen- 
value problem using S(A,B,H) can be solved to find the transmission 
zeroes. The p branches that remain finite approach the LH? transmission 
zeroes, or the mirror image about the imaginary axis of the RHP trans- 
mission zeroes. The zero directions are the vectors associated w.th 
the LHP transmission zeroes. The zero directions associated with the 
mirror image of the RHP transmission zeroes cannot be found using 
this n+m dimension problem. 


! 
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The angles and sensitivities of the optimal root locus can be 
found by applying Theorem 1 and Lcamaa 2 to the namiltonian system. The 
results eore the following: 


Theorem 2. The zmgles on tl»e optimal root locus« for 0 < p < « 
and for distinct are found by 


arg (ds^) 


arg (ds^) 


arg 


arg 




‘o br"^b'^‘ 








H 

Vi 


o 

t 

==i 



0 < p < “ 


H 

Vi 


i » 1/ ...» p 



(25) 

0 < p < ® 

1 « l,...,p 
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Remark. The emgles of departure are found using (25) with 
p-«, and the angles of approach by using (26) with p»0. For p>0, 
p " n> ^md for p« 0, 0 < p < n-m. 


Lemma 3. The sensitivities of distinct closed loop eigenvalues 
to changes in p, for 0 ^ p ^ are found by 


- 1 

1 .H 

r -1 T-< 

0 BR B 


dp 2 

H '^i 


^i 

P 

Vi 

0 0 
u 



0 < p < ® 

1 * 1#...# p 


(27) 


dSi 



dp 

m 

w“z. 


0 < p < ® 


■>» 


i “ 1>...( p 


( 28 ) 
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Rcpark. The conputations for (25**2B) can be reduced by usin^ 
the following identities. First, frccn (23) and (24), it can bo shown 
that 


Second, froa I 8 J , let be the RHP mirror image about the imaginary 
i»xis of s^, and let •» (x^, ^)** be the right eigenvector associated 
with *s^. Then the left eigenvector associated with s^is 


H 


<< 


Xi) 
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I 


) 
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IV. Asymptotically Infinite Properties of the Optimal Root Locus 
In this section we continue to analyse the optimal root locus. 

We review what is known about the asymptotic behavior and then present 
an algorithm which c^m be used to predict the asymptotically infinite 
behavior. 

Review of Known Asymptotic Behavior 

As p 0 the number of asymptotically finite branches is p, where 
0 p ^ n-ci. These branches approach the LHP transmission zeroes of 
the Hamiltonian system. The associated eigenvectors approach zero 
directions, which are part of the zero directions of the Hamiltonian 
system, as explained in section III. The remaining n-p branches group in 
m Butterwoxrth patterns and approach infinity. Let the order of the ith 
pattern by n^. Each of the n^ eigenvalues in this pattern lies on one 

4 

of n^ asymptotes with a distance from the origin approximately equal 
to 



There are n^ right eigenvectors associated with the pattern. These 
span the same subspace of B^spanned by 


4 

A first order Butterworth pattern has one asymptote which coincides 
with the negative real axis. A second order Butterworth pattern has 
two asymptotes which have angles of +45® with the negative real axis. 
In general, an ith order Butterworth pattern has i asymptotes, each 
of which starts at the origin and goes through the LHP solutions s 
of + H)i » 0. 


1 
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00 ^ 

Bv" , ABV*,..., A ^ BV^ . 

The ordering o£ these vectors can be conveniently summarized in terms of 
a multi-index defined in the following way: 

Y « (01, 11, 21,..., 02, 12, 22,..., (n^-lU, ..., 

(n -11m) . 
xa 

Each component ij of Y describes the vector associated witn the m 

Buttcrvorth patterns. 

One special case of the above asimptotically infinite structure 
deserves special notice. When Rank (B*^QB) * m then there are (n-m) 
finite modes find the remaining m infinite modes all form first order 
Butterworth patterns. This is called the "generic" case. For an 
explanation of the word "generic," see (191. The n asymptotically 
infinite eigenvalues lie on the negative real axis an approximate distance 

03 

the origin. Their associated eigenvectors approach Bv^ 
and, hence, the multi-index Y is 


Y * (01# 02# »••# to) 


The asymptotically infinite cigenstructure (s^# i»l#...# mi 
of a generic problem can be readily computed by solving the following 
m-dlmensional eigenvalue problem: 


e» 9 -IT “ 

I (s^ ) I - R B QBJV^ “ 0 


r * l#«.«r 


The resulting solutions will be grouped together in the following way: 


diag • • • f 


OO 03 ^ 

N - V^l 


# ♦ • • # 
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Each of the V* vectors Is a direction and is therefore only unique to 
%dthin a scalar multiple. 

In contrast to the generic problem, the nongeneric case does not 
yield to a similarly simple calculation of its infinite asymptotic 
eigenstructure. For this case, it is necessary to evaluate vectors 

CO QQ 

v^, scalars s^, and also the Buttcrvorth dimensions n^, i - l, n. 

An algorithm for this purpose is provided below. 

^ Algorithn for the Non-Cenerlc Case 

Ondcr our earlier assumptions, the v“ vectors form a basis for 
» . The algorithm presented here decomposes 3R® into these basis vectors. 
This is done in two steps. The first is to compute basis vectors for a 
sequence of subspaces of (defined below) . The second step is to 
use a series of eigenvalue problems to further break down the sub- 
spaces into the v" basis vectors. These same eigenvalue problems compute 
the s**s. 

Let k <_ n-nri-1 be the highest order Butterworth pattern. Define 
the matrices 

J, ■ i=l. k 

and define the subspaces of 

#,0 n 

u w ]R 

» u° n Her 

k * 


‘'k ‘ ‘^k-l 
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Those subspaccs are shovm pictorially in Figure 3* They are nested 
such that 

0 - U^C ... C Uq - ^ , 

and their dimensions satisfy 

- dia - din i » 1,..., k 

k 

B « B . 

i-1 ^ 

A basis for each of the subspaces can be recursively conputed.^ The 
recursion stops at the k step when - 0. Define 

i ■ 0,1, . .., k-1 

to be Batrices whose colunns fona a basis for the subspaces. The 
basis vectors are not imique, and without loss of generality lot 

• 

Though it is not obvious at this point, we note that the number 
of ith order Butterworth patterns is n^. If there are no ith order 
Butterworth patterns then « 0 and - U^. The dimensions of 

are where 



m 


i+l 


m • 


r 


emd is the number of Butterworth patterns of order greater than i* 
When the matrices are computed we have enoiigh information to form 


Numerically this nu y be difficult to do. The decomposition requires deter- 
mination of the kernal of a matrix and the intersection of subspaces. For 
both of the singular value decomposition is a good tool to use [14]. 




Y , the multi-index which lists the orders ol the m Butterworth patterns. 
In the generic case k*» 1, » m, and = 0. 

The next step in the algorithm is to use the 0^ matrices to 
compute N and S . We decompose n" and s” into 

00 00 00 
N - Nj^] 

00 00 00 
S - diag(S^,..., Sj^) . 

00 QJJ 

is an rnxia^ matrix whose colU 2 nns are the •s associated with ith 

order Butterworth patterns, is an m^xm^ diagonal matrix whose 

00 

diagonal elements are the s^'s associated with ith order Butterworth 
patterns. Ke noce that 


= ^i+1 +•••+ 


i“ 0, ..., k-1 


and that in general 


ImN^ Pi 0 for ij^j 

In words, these two equations tell us that one basis for the subspace 

f I 00 

consists of the vectors associated with Butterworth patterns of 
order greater than i, and that in general the ^>*'3 are not orthogonal. 
We define two more sets of matrices: 


q. . 


i If.../ )c 


’'i ■ ‘ ■ 1 " • 


When the matrices are known then the matrices can be computed in 

a straightforward manner. The dimensions of T are I x I 

1 i-1 i-1' 


The 
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sl,.Uic»c o, t.e T, ^ 

.l«..ax„. ,,, „- 

Is made clear In the following Theorea. 


3 _ . The Jordan canonical form of t^ is 


" f*^il «i2l 


A^ 0 


0 0 


!["ii “ix] 


-1 


^ * l#**#f }c« 


(29) 


Aj U a diagoi^l „trl« with poaitlva „al .igaavalaaa, aM 


"i ■ "i-l“ll 


(S.)2 


(30) 


(31) 


Proof* Appendix A. 

Note that in the decomposition ( 29 ) w ta*., 

i«5icion U9) , has as many columns as 

thara ara 1th ordar Batta^tth pattarta. and has aa 

aa thara ata Battat«,tth pattama ,f graataa: than ith padar. a, tha 

dhtanaiana al and raapaativalp, ata and 1^ ^«1 . 

M thara ara no ith ordar hottarvarth pattama than o ani „ u 
not present. ^ 

»ota lurthar that thara ara no raatrictlona on tha a-ltlPUclty 
of tha a", a (thaaa ara poaltiva raal ntrhara and not oloaad loop 
n aluaa). if a^ . tut they ara aasooiated with Buttertorth 
pattama of different ordara than they ara aolutiona to diffarant 
alpanvalua prohlana,!,,. Co„aa,ua„tly thara la no anh.pulty „ m. 
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asscciatod v“ and vectors. However, if s* - s“ and they are associated 
with Buttorworth patterns of the sarae order then wo can only say that 
v“ and Vj form a nonunique basis for a two dimensional subspaco of IR . 
This is known from properties of the Jordan canonical form of T^^ (29) 
and from (30) • 


Example 2 

The algorithm described above is illustrated with the following 
A, B, C - and R matrices i 


0 1 


A ■ 


-5 


0.1 

0 


-4 

0 


‘0 


0 


0.1 


-1 


0 



0 

0 

0 

1 


•65 0 

100 10 



’.01 

I 

• 

O 

R • 

-.01 

.1211 


J 


The asymptotically finite properties are found by solving a generalized 
eigenvalue problem using the system S(ArB#H). The results are 


liol 


1 

•10 
0 

L 6S J 


The as^Tnptotially infinite properties are found by the algorithm 
of Theorem 3. First we find the U subspacos and their matrices U^: 


A 



‘ 30 - 


KB 



KAB 


■[: :] 


since m Xet 




Slnce ker - 0, ^ “ 0* 


The number of first and second order Butterworth patterns are 


« dim - dim ■ 1 

®2 ■ ^2 " ^ • 

The matrices and their Jordan canonical forms are 

’■•■■‘•■■c :]■[: :h: :ic ;] 


-1 


•■■■■■•[:] 


- [ 3 ) 


Tj - {oJrUj^)“^U^(.; 2U^ » (lOOJ 


S - ^«2: - 


a 


s" - aY' ■ ‘iw 


Therefore 


1 1 
1 0 


3 0 

0 10 


Y - (01, 02, 12) 


Each of N represents a direction and is therefore only unique 

to wiih 2 ‘» a scalar multiple. 
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V. Conclusions 

Both the eigenvalue and generalised eigenvalue problems can be used 
to compute angles and sensitivities of multivariable root loci and optimal 
root loci. The generalized eigenvalue problem is superior to use for 
computing the angles of approach. 

The elementary matrices A, B, Q, and R can be used to determine the 

asymptotic behavior of the optimal root locus and the associated eigenvectors. 

A generalized eigenvalue problem can be used to compute S° and X®, the 

asymptotically finite properties. A subspace decomposition of the control 

space and a scries of eigenvalue problems can be used to compute 
00 00 

S , M , and Y; the asymptotically Infinite properties. 

We are hopeful that a similar type of subspace decomposition can be 
used to determine the asympotically infinite behavior of arbitrary 
multivariable root loci. The extension of the present method is diffi- 
cult, however, because we do not in general have the synmetry of closed 
loop eigenvalues about the Imaginary axis forced by the opttoal Hamiltonian 


system. 
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Appendix A - Proof of Theora 3 

The proof is by Induction and uses the fact that all closed loop 
eigenvalues s^ and closed loop vectors sust satisfy 

tPR + «’’(-9jl«tSj)IV^ - 0 J . 1 

vrfiere 

♦(s) « H(sX-A)“^B 


H(il + ^ A +...)B 


1-1 s ^ 


Equation (A.l) is derived in [20]. It can also be derived by 
manipulations of (23). An expanded version of ^'^(-s)'Ks) is 


♦^(-s)«(s) - £ ~ ^ (-l)^j/j, 1 • 

i-2 Ls ^ 


The first step in the induction proof is to show that the theorem 
is valid for and We assume without loss of generality that 
first order Buttervrorth patterns exist. Equation (A.l) can be 
rewritten 


PR + ~ b'^QB + 
"j 



j — li...f n . 


(A. 2) 
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A8 the tern doainates (for the asynptotically infinite eigen- 


v&XuQ3 s^) and (A«2} can bo roi^ribben 


a — 1 T 

jl - R B QB)V^ " 0 


i ** B 


X. «» psT 




CA.3) 

(A.4) 


The eigenvalues of r"^b^QB are real and nonnegative. (This is be- 
cause the eigenvalues are the sane as those of 

a natrix of the form x\, which is known to have real and nonnegative 
eigenvalues), whan > o we can use (A.4) to solve for s^. The I.HP 

solution is s^ « which is a first order Butterworth pattern 

with - xJ/2. Therefore s" as given in Theorem 3 is valid. h“ is also 
valid, be’csuoe from (A. 3) we see that the vectors associated with 
first order Button«rth patterns are eigenvectors of r”^b"^QB. 

The Vj vectors associated with Butterworth patterns of order greater 
than one form a basis for the kernel of r'^b^qb, which is tf^. Heuristically 
speaking, these v" vectors are not “trapped- by the s~^ term of (A. 2). 

The next step in the induction proof is to assume that s“ and 

00 X— X 

“i-l Chat s“ and m" are valid, tf s“_^ and m“ ^ 

are valid then the vectors associated with Butterworth patterns of 
order ^ i form a basis for Therefore for each of these vectors 

there exists an vector such that v“ - . Substitute this into 

(A.l) to get 


[pR + 0 




i-I 



Multiply on the left to got 




O^.ltpR ♦ - 0 j » 1,.... 1 


i-l‘ 


After some algebra this reduces to 


* t-i)‘ rh-’I-AVi * ^ "j - 0 • 


-2i 


i Vi 

As the term dominates and we get 

i * f 


(\^l - T^)01^ - 0 


i-1 


(A. 5) 




(A. 6} 


The eigenvalues X^ are real and nonnegative, for the saste reasons as 
for the i«l case* When > 0 we can solve for q^, and the LUP solutions 

CO 1/2 

are recognized as an ith order Butterworth pattern with s^ * X^ . 

00 

Therefore of Theorem 3 is valid. The eigenvectors 6}^ of T^ associated 
with the nonzero eigenvalues X^ are the columns of and therefore 

CD 00 

» ^i-l'^il* associated with Butterworth patterns of 

-2i 

order greater than i axe not **trapped** by the term, and therefore 

they lie in This completes the proof. 
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